Phase Transition in a Traffic Model with Passing 
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We investigate a traSic model in which cars either move freely with quenched intrinsic velocities 
or belong to clusters formed behind slower cars. In each cluster, the next-to-leading car is allowed to 
pass and resume free motion. The model undergoes a phase transition from a disordered phase for 
the high passing rate to a jammed phase for the low rate. In the disordered phase, the cluster size 
distribution decays exponentially in the large size limit. In the jammed phase, the distribution of 
finite clusters is independent on the passing rate, but it accounts only for a fraction of all cars; the 
"excessive" cars form an infinite cluster moving with the smallest velocity. Mean-field equations, 
describing the model in the framework of Maxwell approximation, correctly predict the existence of 
phase transition and adequately describe the disordered phase; properties of the jammed phase are 
studied numerically. 



PACS numbers: 02.50-r, 05.40. +j, 89.40-|-k, 05.20.Dd 



I. INTRODUCTION 

TrafRc flows on single-lane roads with no passing ex- 
hibit clustering since queues of fast cars accumulate be- 
hind slow cars. These clusters form and grow even when 
car density is small. The initial analysis of cluster forma- 
tion was carried out in the earlier days of traffic theory 
0, and this subject continued growing ever then [||-||. If 
passing is introduced, the clusters may stop gro wing af- 
ter reaching a certain size. Indeed, previous work |l0f-|l|] 
indicated that after a transient regime a steady state is 
reached. The models of Refs. |10|-[l2[ assume that any 
car in a cluster can pass the leading car and the passing 
rate is independent on the location of the car within the 
cluster. This is certainly an oversimplification of the ev- 
eryday traffic scenarios. The complementary case when 
only the next-to-leading cars can pass is also an ideal- 
ization, yet it is closer to reality. Below we show that 
the latter model is also richer phenomenologically as it 
undergoes a dynamical phase transition. 

We first comment on possible theoretical approaches. 
A mean-field theory is the primary candidate, and we be- 
lieve that it may be very good, perhaps even exact, since 
clustering and passing mix positions and velocities of the 
cars. The Boltzmann equation approach is an appropri- 
ate mean-field scheme, and in our earlier work pO|JTl[ | we 
indeed used it. However, the present model, where only 
the next-to-leading car is allowed to pass, is significantly 
more difficult than the model where passing was 

possible for all cars. Indeed, it appears impossible even 
to write down closed Boltzmann equations for the dis- 
tribution functions like P{v,t) and Pm{v,t), the density 
of all clusters moving with velocity v, and the density of 
clusters of m cars, respectively. Therefore our theoreti- 
cal analysis is performed in the framework of the Maxwell 
approach. This scheme simplifies "collision" terms by re- 
placing the actual collision rates, which are proportional 
to velocity difference of collision partners, by constants. 
Despite this essentially uncontrolled approximation, the 
Maxwell approximation is very popular in kinetic theory 



|i3[ and it has already been used in traffic [|2[. 

The important feature of our model is quenched disor- 
der, which manifests itself in the random assignment of 
intrinsic velocities. Road conditions (construction zones, 
turns, hills, etc.) present another source of quenched 
randomness in real driving situations p^ , which is ig- 
nored in our model. Quenched disorder significantly af- 
fects characteristics of many- particle systems, especially 
in low spatial dimensions |15| . This general conclusion 
applies to the present one-dimensional traffic model as 
we shall show below. 



II. MAXWELL APPROXIMATION 

We now formally define the model. Free cars move 
with quenched intrinsic velocities randomly assigned from 
some distribution Po(v). When a car or a cluster encoun- 
ters a slower one, it assumes its velocity and a larger 
cluster is formed. In every cluster, the next-to-leading 
car is allowed to pass and resume driving with its in- 
trinsic velocity. The rate of passing is assumed to be a 
constant. Thus clusters move and aggregate determinis- 
tically, while passing is stochastic. The system is initial- 
ized by randomly placing single cars and assigning them 
uncorrelated intrinsic velocities. 

Within the Maxwell approach, the joint size-velocity 
distribution function (the density of clusters of size m 
moving with velocity v) Pm{v,t) obeys 



dPm{v,t) 

dt 



= 7(1 - 5m,l)[Pm+l{v,t) - Pm{v,t)] 

+ l5raA[N{v,t) + P2{v,t)] ~ c{t)P^^{v,t) 
nOO 

+ dv' mv',t)P,{v,t). (1) 

^ V - I ■ 

i+j—m 



Here 7 is the passing rate, so terms proportional to 7 
account for escape, while the rest describes clustering. 
The escape terms are the same within Boltzmann and 
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Maxwell approaches, and they are actually exact. The 
collision terms are mean-field by nature, and they are 
different in the Boltzmann and Maxwell approaches. For 
instance, in the Boltzmann case, the integral term must 
involve v' — v. Eqs. (|f|) also contain c(t), the total cluster 
density 



c{t) 



E 

3>1 



dv Pj {v,t), 



(2) 



and N{v,t), the density of clusters in which the next-to- 
leading car has intrinsic velocity v. This N{v,t) causes 
the major trouble since it cannot be expressed through 
Pj{v,t). One might try to close Eqs. (|l]) by introducing 
Fk{v,v',t), the density of clusters moving with the ve- 
locity v' whose k^^ car has intrinsic velocity v. Clearly, 
N{v,t) = dv' F2{v,v' ,t), and it appears that equa- 
tions for Fk{v, v' , t) are closed. A more careful look, how- 
ever, reveals that the governing equation for F2(v,v',t) 
includes three-velocity correlators. 

Thus, at the first sight, the Boltzmann and Maxwell 
approaches appear to be equally incapable of providing 
closed equations for the joint size-velocity distribution 
function. Still, the Maxwell framework has an advantage 
that it does provide a closed description on the level of 
the cluster size distribution. Indeed, integrating Eqs. (|l|) 
over velocity and defining Pm{t) = dv Pm{v,t), we 
find that the cluster size distribution Pm{t) obeys 



dP,n 

dt 



i+j—m 



for m > 2, and 
dPi 

~dr 



l[P2- Pl+C\~cPi. 



(3) 



(4) 



Besides this formal derivation of Eqs. @-(||) by direct 
integration of Eqs. (|l]), it is possible to obtain these equa- 
tions by enumerating all possible ways in which clusters 
evolve. For instance, consider Eq. (^. Collisions re- 
duce the density of single cars, and the collision rate is 
clearly equal to c{t), as it is velocity- independent in the 
framework of the Maxwell approach. The escape term in 
Eq. (^) is understood by observing that the rate of return 
of single cars into the system is equal to 



J>3 



= l[P2- Pi 



Here P2{t) is singled out since passing transforms it into 
two single cars while an escape from larger clusters pro- 
duces only one freely moving car. 

Eqs. are closed. Mathematically similar equa- 

tions were investigated previously in the context of the 
aggregation-fragmentation model Therefore, we 

merely present essential steps of the analysis. Restricting 



ourselves to the steady state and introducing notations 
Pni = iFm, Coo = 7^: ^6 rccast Eqs. into 

FF„,^F„,+^-F,r, + 5„^sF+\ E ^-^^y 

i+j—7n 

These equations should be solved together with the con- 
straints Y.ra>l P^ri = Cqo and Y.ra>l ^ 1' i'^-; 



(6) 



Note that the sum X)m>i "^^'"('•) obviously constant 
due to car conservation. The constant is equal to the 
initial concentration cq as cars were initially unclustered. 
Here and below we always choose cq = 1. 

As in Ref . pql , we introduce the generating function 



T[z)^\^{z"^-l)F„ 

m>l 

This generating function obeys 



1 



1- z 



2 z 
with the solution 

z - 1 



■J" 



[l-zf 



P = 0, 



T{z) = ^— |l - Vl - 2zf| 



(7) 



(8) 



(9) 



The steady state solution exists only when the gen- 
erating function is real for all the Q < z < 1. Hence, we 
require that 2F < 1. Assuming that this condition is sat- 
isfied, we expand the generating function in the powers 
of z to obtain the steady state concentrations: 

(2^)™ f r (m - i) r(m-|-i)l , ^ 
p^ = \ i <; U- ~ 2F U- \ ■ (10) 



20F I r(m-l-l) 



r{m + 2) 



This solution is still incomplete as we have not yet de- 
termined F. To find F we use the sum rules (H). The 
first sum rule is manifestly obeyed, while the second sum 
rule yields J^mFm = dT/dz\^=x = 1 - y/l~2F = 7-1. 
Thus, F = -^r-, which translates into Cqo = 1 — I/27. 

The steady state solution ( p^ ) exists for sufficiently 
high passing rates, 7 > 7c = 1. For 7 > 1 and large m, 
the steady state size distribution simplifies to 



Pn 



(11) 



with C = (47r) ^/"^^ ^(7 — 1)^. Apart from a power-law 
prefactor, the size distribution exhibits an exponential 



decay, P„ 



, in the large size limit. The charac- 



teristic size diverges, m* ~ (7 — 1)" as the passing rate 
approaches the critical value 7c 1. In the critical case, 
the size distribution has a power-law form 



p — 



3 r (to 



-5/2 



4v^ T{m + 2) 



(12) 
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Let now the passing rate drops below the critical value 
(7 < 7c)- Since F cannot grow beyond Fc = 1/2, it stays 
constant. Therefore, Fm is given by the same Eq. jl^ ) 
as in the critical case, and the cluster size distribution 
reads Pm = "fFm- This implies Coo = 7/2, i.e., the sum 
rule ^ Prn ~ Coo is valid. The second sum rule is for- 
mally violated: ^mPm = 7 7^ 1, i.e. the cluster size 
distribution (jlj) accounts only for the fraction of all the 
cars present in the system. The only possible explana- 
tion is the formation of an infinite cluster that contains 
all the excessive cars. The second sum rule then shows 
that 1 — 7 of all the cars in the system are in this infinite 
cluster. 

Thus within the framework of the Maxwell approxima- 
tion, our traffic model displays a phase transition which 
separates the disordered and jammed phases. The steady 
state cluster concentration has different dependence on 
the passing rate for these two phases: 



_/ 1-1/27, 
"17/2, 



7 > 1; 
7 < 1. 



(13) 



In the disordered phase, the size distribution decays ex- 
ponentially in the large size limit. In the jammed phase, 
Pm has a power law tail and in addition there is an infi- 
nite cluster which contains the following fraction of cars: 



0, 



7, 



7 > 1; 
7 < 1. 



(14) 



This phase transition is similar to phase transitions in 
driven diffusive systems without passing [^|-0J^ and to 
phase transitions in aggregation-fragmentation models 
p6|-p^ . Also, the mechanism of the formation of the infi- 
nite cluster has a strong formal analogy to Bose-Einstein 
condensation Pp7|. 

Turning back to the joint size- velocity distribution (|^), 
we note that the lack of an exact expression for N{v) in 
terms of Pm{v) does not mean the lack of a mean-field re- 
lation between these quantities. Indeed, the density N{v) 
of clusters in which the next-to-leading car has intrinsic 
velocity v, can be written as 



N{v) 



f dv'Y^p.w 

■Jo ,>2 



Civ) 



Cdv"C{v") 



(15) 



Here '^^^2^3^^') density of "true" clusters (i.e., 

freely moving cars are excluded) moving with velocity v' . 
Then, C{v) = Pq{v) — P{v) is the density of cars with 
intrinsic velocity v which are currently slowed down, i.e., 
they are neither single cars, nor cluster leaders. Assum- 
ing that the velocities of cars inside clusters are perfectly 
mixed, C{v)/ dv" C{v") gives the probabihty density 
that the next-to-leading car in a true ti'-cluster has the 
velocity v. The product form of Eq. ( |l5| ) reveals its mean- 
field nature, which is consistent with the spirit of our the- 
oretical approach. One can verify that Eq. ([l^) agrees 
with the sum rule / dvN{v) = J2j>2^j^ thus providing 
a useful check of self-consistency. 



Although Eqs. (|l|) with N{v) given by (y_5|) seem very 
complex even in the steady-state regime, several conclu- 
sions can be derived without getting their complete solu- 
tion. We first simplify Eqs. (|l|) by introducing auxiliary 
functions 



(v) = / dv'Pmiv'). 



(16) 



By inserting P,„ — — into the Eqs. (|i|), integrating 
resulting equations over v, and using the boundary con- 
ditions Qm{v = 00) — 0, we find 

7 [<3m+l(w) - Qm(w)] - cQm(w) + i ^ Qi{v)Qj{v) 

= S,ni7q{v), (17) 



with 



qiv) = -Qiiv) ~ \ dv'Niv') 



(18) 



Eqs. ( p^ ) are almost identical to the Eqs. (^)-(^, the 
velocity is just a parameter. Consequently, we anticipate 
qualitatively similar results, Qm{v) ^ ^-3/2g-m/m ^ 
and 



(19) 



with the characteristic size m*(v,j) dependent on both 
velocity and passing rate. Our more rigorous generat- 
ing function analysis, performed along the lines described 
above, confirms the asymptotic form (Eoh. 



III. SIMULATIONS 

Now let us examine what conclusions obtained within 
the Maxwell approach are relevant for the original model. 
We first re-derive the condition for the phase transition in 
the complete velocity-dependent form. Let us consider a 
system of reference with the origin moving with the slow- 
est car. We assume that the system is sufficiently large 
for the slowest car to have negligible velocity. We com- 
pare the total flux of cars clustering behind this slowest 
car, '^mPm{v)m, to the rate of escape, 7. Here {v)m 
is an average velocity of a cluster of size to. When the 
rate of escape becomes less than the rate of accumulation 
of the cars, the cluster behind the slowest car (analog of 
the "infinite cluster" for finite systems) grows to remove 
the excessive cars from the system. Hence, the phase 
transition point 7c is defined as 



E 

m>l 



■mP„ 



Ic- 



(20) 



For the Maxwell model, where — 1 for all to, 

Eq. ( pO| ) reduces to ^ mPm. = 7c = 1 as obtained above. 
Since large clusters usually form behind slow cars, {v)m 
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is a decreasing function of the cluster size m. In particu- 
lar, {v)m is always smaller than the average car velocity 
(v), implying > < 1. 

For a rough estimate of {v)m, consider a cluster of m 
cars and assume that intrinsic velocities of the cars in the 
cluster are independent. The leading car has the minimal 
velocity, so the size-velocity distribution reads 



Pm{v) ~ mPo{v) 



dv'Poiv') 



Pr, 



(21) 



7O 



□ r 



v-0 

V 



o 



'V 



^0 



For concreteness, let us consider intrinsic velocity distri- 
butions which behave algebraically near the lower cutoff, 
Pq{v) ^ yf^ as w — > 0. Then for large clusters we get 



Pm{v) ^ Pm exp (- 



)■ 



(22) 



This implies that the average cluster velocity {v)rn scales 
with m according to {v)m m~'^/'^^^^^\ and hence 
7c ~ ^ m^/^'^+^-'Prn- We conclude that the phase tran- 
sition does exist in the original model, although its lo- 
cation is shifted towards lower passing rate compared to 
the Maxwell model prediction. This shift is especially 
significant for small n (fi > —1 from the normalization 
requirement). 

To check the relevance of other predictions of the 
Maxwell approach, we performed molecular dynamics 
simulations. We place N = 20000 single cars onto the 
ring of length L = N, so that the average car den- 
sity is equal to one. Initial positions and velocities 
of cars were assigned randomly. We considered linear 
P„{v) = |w (0 < w < 3/2), exponential Po{v) = er" , 
and Po{v) = (27rti)~^/^e^'"/^ velocity distributions, which 
correspond to fi — 1, 0, —1/2 for the small-w asymptotics. 
All these three distributions have the average velocity 
equal to one. 

In Fig. 1, we plot ln[m^/^Pm] vs. m for the above three 
velocity distributions. We take 7=1 which, as we con- 
cluded before, lies above the phase transition point jc- 
We expect the system to be in the disordered phase with 
P,n being expressed by Eq. (|ll]). For the exponential 
and Poiv) — {2T:v)^^^^e~'"^^ intrinsic velocity distribu- 
tions, there is a good agreement with the prediction of 
the Maxwell model (11); for the linear velocity distribu- 
tion, there are some deviations for small m, but for large 
m the agreement is satisfactory. The slopes of the plots 
decrease with fi. Taking into account that at the point of 
the phase transition the slope equals to zero, this quali- 
tatively confirms that 7c gets smaller when ^ decreases. 



Fig. 1. Plot of in [m^'^^Pm] vs. cluster size m in the high 
passing rate regime (7 — 1) for linear (□), exponential (o), 
and Po(v) = (27ri;)~^''^e~"''^ (V) initial velocity distributions. 

Plots of Pm vs. m for intrinsic velocity distributions 
Po{v) = e-" and Po{v) = (27rw)-i/2e-^/2, with passing 
rate 7 = 0.005 well below the phase transition point, are 
shown in Fig. 2a and Fig. 2b, respectively. The cluster 
size distribution clearly consists of two regions: almost 
power-law tail for smaller m and several separate peaks 
for larger m. These peaks correspond to the fluctuating 
size of the infinite cluster, while the power-law tail de- 
scribes the regular part of Pm- The apparent exponent 
r of the power-law region Pm ~ m^'^ slightly varies for 
different passing rates and Po{v), though it remains con- 
fined between 3/2 and 2. It is definitely different form 
the value 5/2, predicted by the Maxwell model (p^). The 
measured values of r would make the total amount of cars 
in the system divergent, ^ mPm — > 00, so the power-law 
region ends with an exponential cutoff at large m. 

We now comment on the relationship of our model to 
earlier work. On the mean-field level, our model is similar 
to the models of Refs. Jl6| , p7t . On the level of the pro- 
cess, our model reminds an asymmetric conserved-mass 
aggregation model (ASCMAM) |17| where clusters un- 
dergo asymmetric diffusion, aggregation upon contact, 
and chipping (single-particle dissociation). Of course, 
our model is continuum while the ASCMAM is the lat- 
tice model. More substantial difference between the two 
models lies in the nature of randomness - in our model 
intrinsic velocities are quenched random variables, while 
in the ASCMAM dynamics is the only source of random- 
ness. Nevertheless, the phenomenology of the two mod- 
els appears to be quite similar. In particular, the ASC- 
MAM undergoes a phase transition, and in the jammed 
phase, the cluster size distribution exibits a power law 
decay with the exponent close to 2 [0 . We should stress 
that in the jammed phase, we have not reached a scale- 
free critical state which must have the exponent r > 2. 
Maybe quenched randomness does not allow the system 
to organize itself into a truly normalizable critical state. 
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Other possible explanation relies on large fluctuations in 
disordered systems, i.e., our system was not large enough 
to ensure self- aver aging. 




Fig. 2a. Plot of the steady state cluster size distribution 
Pm in the low passing rate regime (7 = 0.005) for the expo- 
nential initial velocity distribution. 



the jammed phase, where the cluster size distribution be- 
comes independent on 7 and the infinite cluster is formed. 
Within the framework of Maxwell approach, which plays 
the role of the mean-field theory in the present con- 
text, we have shown that the model admits an analytical 
solution. We have argued that the Maxwell approach 
correctly predicts the existence of the phase transition 
and adequately describes the properties of the disordered 
phase which arises when the passing rate is high. For 
the jammed phase, the Maxwell approach correctly pre- 
dicts that the system stores excessive cars in the infinite 
cluster and organizes itself into some kind of a critical 
state. However, the Maxwell approach cannot quantita- 
tively describe other properties of the jammed phase. It 
would be interesting to design a more accurate theoretical 
approach which would allow to probe the characteristics 
of the low passing rate regime analytically. Some prop- 
erties of the jammed state appear similar to the proper- 
ties of the jammed state of a lattice model of Ref. jl^] 
which includes an asymmetric lattice diffusion, aggrega- 
tion, and fragmentation. It would be interesting to gain 
a deeper understanding of the relationship between these 
models, and whether the quenched disorder is the main 
source of difference. 



We are thankful to E. Ben-Naim and S. Redner for dis- 
cussions, and to NSF and ARO for support of this work. 




Fig. 2b. Plot of the steady state cluster size distribu- 
tion Pm in the low passing rate regime (7 = 0.005) for the 
Poiv) = (2Tvv)~^^^e~^^^ initial velocity distribution. 



IV. CONCLUSION 

In this paper, we have investigated the model of traf- 
fic that involves clustering and passing of the next-to- 
leading car. Despite the fact that it is one of the sim- 
plest (if not the simplest) possible continuous model of 
one-lane traffic with passing, the model has rich kinetic 
behavior. Depending on the passing rate 7 the system 
organizes itself either into disordered phase where den- 
sity of large clusters is exponentially suppressed, or into 
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